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Abstract. We study a single-flip dynamics for the monotone surface in (2 + 1) di- 
mensions obtained from a boxed plane partition. The surface is analyzed as a system 
of non-intersecting simple paths. When the flips have a non-zero bias we prove that 
there is a positive spectral gap uniformly in the boundary conditions and in the size 
of the system. Under the same assumptions, for a system of size M, the mixing time 
is shown to be of order M up to logarithmic corrections. 
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1. Introduction, model and results 

Consider a surface in 2 -|- 1 dimensions defined by non- negative integer heights ix,y^ 
where x,y £ such that £x,y ^ ^x+i,y and ix,y ^ ^x,y+i for all x,y G When 
ix,y = for all but finitely many x, y this is called a plane partition, the two-dimensional 
generalization of an ordinary partition (Young diagram) . When the surface is such that 
^x,y ^ c and ix,y = when either x ^ a or y ^ 6, for some integers a, b, c, then it defines 
a boxed plane partition, or a plane partition in the box a x b x c. As we shall see, a 
convenient representation of a boxed plane partition is obtained by considering a system 
of non-intersecting simple lattice paths. Other well known equivalent characterizations 
are the perfect matchings (or dimers) configurations on a subgraph of the honeycomb 
lattice and the lozenge tilings of an hexagon. 

A continuous time flip dynamics of a plane partition in the box a x b x c is defined 
as follows: Every point (x, y) in the rectangle Ra,b = {0, . . . , a — 1} x {0, . . . , 6 — 1} 
is equipped with an independent, rate 1, Poisson clock. When (x,y) rings we flip an 
independent {0, 1} coin X; if X = 1 we replace £x,y by i'^^y = lx,y + 1 if allowed; if 
X = we replace lx,y by I'^ y = ix,y — 1 if allowed. When the coin is unbiased the flip 
dynamics converges to the uniform distribution over plane partitions in the box axbxc. 
The latter is known to exhibit non trivial limiting shape or "arctic circle" phenomena 
in the limit of large parameters a, b, c, cf. [3], [7] and references therein. We refer to [2] 
for the connection with the Wulff crystal in the low-temperature 3D Ising model. It 
is an interesting open problem to determine the speed of convergence to the uniform 
equilibrium measure. We refer to ^ and references therein for the various polynomial 
bounds known so far and for the conjectured diffusive behavior of the spectral gap. 

In this paper we shall analyze the biased case. When X = 1 with probability p ^ 
then the flip dynamics converges to the probability measure fia on plane partitions in 
the box axbxc such that every configuration £ has a weight proportional to e-2aVoi(£) ^ 
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where Vol(^) stands for the volume y ^x,?/ under the surface, and e = — p). 
By symmetry, fia is equivalent to fi-a and we shall restrict to positive values of a (or 
p < 1 /2) . In this model there is no critical value of a and one has a localized surface 
(i.e., Vol(£) = 0(1)) for all a > 0. This follows from the fact that the number of plane 
partitions i such that Yol{£) = v is e*^^"^ ^\ see Section [2] for more details. 

It has been recently shown that a direct coupling argument allows one to prove 
that, if a is sufficiently large, then uniformly in the size of the box one has a positive 
spectral gap and a mixing time of order M = max{a, b, c}, see jB]. Below we shall prove 
that this actually holds for all a > 0, up to a logarithmic correction in the mixing time 
upper bound. As in [B] we use path-coupling arguments with an exponentially weighted 
metric. However, in the case of small a > one of the novelties is that these arguments 
can only be applied to suitable coarse-grained versions of the process. In the proof of 
the spectral gap estimate we compare the single-flip process to two auxiliary coarse- 
grained dynamics. Another novel ingredient is a roughly deterministic description of 
the process at large scales which allows one to obtain the mixing time estimate. 

It is worth pointing out that the measure //q,, a > 0, has a natural extension to 
infinite boxes, i.e. as a measure on plane partitions without any box constraint. Our 
spectral gap estimates imply that for all a > this extended measure has a positive 
spectral gap. 

1.1. The model. We first formulate the model in terms of configurations of non- 
intersecting paths and then describe the mapping needed to obtain boxed plane parti- 
tions. 

Let k, n, h be integers, such that k,n ^ 1 and h £ {—n, . . . , n}. We consider the set 
n collections of polymers described as follows. Each polymer is a one-dimensional 
nearest neighbor path of length n which starts at height and ends at height h, and 
there are k ordered paths: a configuration t] £ Q,^^ is characterized by integer heights 

r]x G Z, j = 1, . . . , /c and x = 0, . . . ,n satisfying the constraints: 

ryj^) = , r/(f ) = h , Vr/(J) (x) := rj^Jl-^ - G {-1, +1} , and 

Vx^>Vx^^K l^j^A;-l,0^2;^n-l. (1.1) 

The set Q,^ ^ is non-empty if h and n have the same parity. Given a > 0, the equilibrium 
measure IJ' = IJ'k n h defined by 



exp ( a }2^=l 22x=o 

Kv) = 

where 



M^) = ^ ^ V&^k,n, (1-2) 



k 



^ = ^M,/.= E exp a^^^r^W 

veni„ \ 1=1 x=o / 

is the normalizing constant. 

For every ^, a G we will write simply ^ ^ a when ^x ^ (^x, ^ x ^ n. For a 

given pair ^ ^ o" we consider the subsets 

^C,<x = {r/ef^t' and,?W^a}. 
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We call and a the ceiling and the floor respectively, for intuitive reasons, see Figure 
[TJ Given a ceiling ^ and a floor o" ^ o") we define the equilibrium measure n^'^ by 
conditioning on £^^,0-: 

//^''^ = /i(-|i?^,.) . (1.3) 
Note that we may write more explicitly 

^C.'^(^) = ! J^VSEiA. (1.4) 



where 



l^j<fcls£x^ n-1 

is the "total volume between ^ and r/", and 



^ e-2"^'(''). (1.6) 




Figure 1. A configuration r/ G -E^.o-. Here A; = 6, n = 18, /i = 6. For 
graphical convenience ry(*) , the i-th component of 7], has been shifted by 
—i units in the vertical direction. The top path is the ceiling ^ while 
the bottom path is the floor a shifted by —{k + 1) units. 

We will denote by A := A^""'^^ the maximal one-polymer configuration in ^i^'- 

_j X for ^ X ^ (n + /i)/2 . . 

'''~\{n + h)-x for (n + /i)/2 < x ^ n. ^ ' 

We will also use the notation V := V^"''^-' for the minimal one-polymer configuration: 

V .= / for ^ X ^ (n - h)/2 

^' \ {h - n) + X for {n - h)/2 !^ X ^ n. ^ ' ' 

Note that {n + h)/2 is an integer since n and h have the same parity, and that V = — A 
if /i = 0. 
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1.1.1. From non-intersecting paths to plane partitions. Suppose first that C = cind 
cr = V so that there is no further constraint on S ^fc„- Then there is a bijection 
between ^ and the set of plane partitions in the box a x b x c with c = k, a + b = n, 
b — a = h (i.e. a = {n — h) /2, b = {n + h)/2). The map is best explained informally as 
follows. From the configuration r] we get the stepped surface i{r]) by adding layers of 
height 1 to the basis rectangle Ra,bj according to the paths r]^^^ chosen in reverse order. 
Namely, r/*^'^-' is the first layer. On top of that we put the second layer jy^^^^) and so on. 
This defines a bijection, see also Section [4.1.11 for more details. 




Figure 2. The stepped surface corresponding to the paths in Figured) 

For any ^ G $7^ j^., let ^ ^ ^fc n denote the configuration such that = ^, i = 1, . . . ,k. 

Also, let £{^) denote the associated stepped surface. Then it is easily seen that, for any 
^ ^ CT, the map described above gives a bijection between the set E^ ,j and the set of 
plane partitions £ in the box a x b x c (as above) such that £{^) ^ i ^ ^(''')> see Figure 
[2j Note that here inequalities are reversed with respect to the polymer representation. 

1.2. Results. The following heat bath dynamics for configurations of lattice paths is 
easily checked to be equivalent to the fiip dynamics discussed in the introduction. We 
define the continuous time Markov chain on the set -E'g,(ji for a given pair ^, a G ^in^ 
^ ^ o", as follows. At each (i,x), i = l,...,k and x = l,...,n — 1, there is an 
independent rate 1 Poisson clock. When (i, x) rings we update the height rjx^ with a 

~(i) 

new height r]x sampled according to the conditional distribution 




where r]^^^ = ^ and r^^'^^^) = a. The Dirichlet form of this process is given by 

k n— 1 

^(/) = EE/^^'''[Var..(/)] , (1.9) 

i=l x=l 

where Vai^i^xif) denotes the function 
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and / : ^ — )• M denotes an arbitrary function. The spectral gap is given by 

£{f) 



gap^o; ) = mm 



/ Var(/)' 

where Var(/) = fJ,^''^{f'^) — /^^''^(/)^, and the minimum ranges over all / : 17^ ^ — )• M 
such that Var(/) 7^ 0. 

Theorem 1.1. For any a > 0, there exists c{a) > such that, uniformly in k,n £ N, 
\h\ ^ n and ^ ^ a G ^i^n- 

gap(a) ^ c(q) . (1.10) 

The estimate in Theorem 11.11 is already known to hold when k = 1 (see |li Th. 4.3] 
or, for an alternative proof, [6]). 

Our second result concerns the mixing time of the Markov chain which, we recall, is 
defined as 

Trnix = inf{t > : max ||Pt(r/, •) - ^||var ^ l/(2e)}, (1.11) 
with II • 1 1 var denoting the total variation norm: 

Ik - Z^'llvar := ^ ^ W{V)-'^'{V)\- (1-12) 

Pt{r]-, •) is the law, at time t, of the Markov chain started from r/ at time zero. 

Theorem 11.11 implies that the mixing time of the Markov chain defined above in the 
case 2A: = n, /i = 0, ^ = A, cr = V is 0{n^). This is a simple consequence of the 
well-known inequality 

Tmix < gap"^ ^1 - log ^mm^(r?)^^ , (1.13) 
(see also Lemma 14.31 below) . We can however prove: 

Theorem 1.2. Let M = max(n,A;). For every a > there exists C{a) < 00 such that 
uniformly in the choice of ceiling ^ and floor a, T^aix ^ C{a) M(log M)^. 

We mention that the conjectured behavior is 0{M), without logarithmic corrections. 
This bound is known to hold if a is sufficiently large [6J . 

For simplicity we have stated these results for a positive constant bias a > 0, but 
there is no difficulty to extend them to the case of a non-homogeneous bias Ox^y on 
each column of the stepped surface ix,y, provided that there exists oq > such that 
Oix,y ^ ao for all x, y. Alternatively, one could place a non-homogeneous bias a(j j.) 
associated to each polymer i and position x in the definition of the measure (|1.2p . 

The rest of the paper is organized as follows. In Section [2] we provide some prelim- 
inaries and prove a couple of key equilibrium estimates to be used in the proof of the 
main theorems. The latter is given in Section [3] (Theorem [TTT]) and Section U] (Theorem 
O). 
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2. Preliminaries 

2.1. Particles and vacancies. Each polymer rj^^^ can be characterized by the posi- 
tions of its positive increments, also called particles. More precisely, for every 1 ^ j ^ k, 
let x'f^ = x'f\r]), i ^ 1 denote the position of the i-th positive increment in the j-th. 
polymer, defined recursively by: 

x^/^ = min{x € {0, ... ,n - 1} : Vr]^^\x) = +1} , . . . 

x^J-^^ = mm{x > x^p : Vr]^^\x) = +1} . 

Note that, given h, the k polymers in the configuration rj £ ft^^ all have the same 

number = N{h,n) = (n + h)/2 of particles (they all have the same length n, the 

same starting point and the same end-point h). Observe that, because of the order 

(i) (7+1) 

constraint, particles obey the following relations: xf ^ , i = 1,...,N and 

j = 1, . . . ,k — 1. We often write x or x(r]) for the collection of particle positions of a 
given configuration rj. 

The set of vacancies for the polymer rj^^^ is defined as the set of points in {0, . . . , n— 1} 
which do not contain particles. Of course, the number of particles plus the number of 
vacancies for r/ G fij^^ equals n. 

2.2. Monotonicity. Trajectories of the Markov chain corresponding to distinct initial 
conditions and/or distinct boundary constraints can be realized on the same probability 
space by a standard coupling argument. This is a straightforward generalization of the 
argument for a single polymer, see O Section 2]. It follows that the Markov chain 
enjoys the following useful monotonicity property. If ri^''^{t; (") denotes the evolution of 
the surface with ceiling ^ and floor a at time t and with starting configuration ( at 
time 0, then almost surely one has 

77«''^(t;C)^/''^'(t;0, (2.1) 
whenever ^ ^ cr ^ a' and C ^ C'. Here for two systems of polymers r], C S „ we 

use the convention that ?y ^ C means r/^*) ^ i^^*) for all i. 

Let E denote expectation with respect to this global coupling P. Using the notation 
IE[/(r/^''^(t;C))] = PtfiC), f ■■ ^k,n then implies that for every fixed t ^ 0, 

the function Pff is increasing whenever / is increasing, where a function / is called 
increasing if /(?]) ^ /(C) ^^r any rj,( such that ^ C- Useful inequalities for the 
equilibrium measures can be derived from this. For instance, taking the limit t — )• 00 
in ()2.ip yields the inequality 

M^''^(/)^/''^'(/), (2.2) 
for any increasing / and any ^ ^ ^',a ^ a'. We will often use one form or another of 
the inequality ()2.2p without explicit reference. 

2.3. Tightness of the excess volume and decay of correlation. Here we prove 
some equilibrium results concerning the exponential decay of spatial correlations, and 
the exponential tightness of the "excess volume" y^(r/). While our main aim is to 
provide the necessary tools for the proof of Theorems 11.11 and II. 2^ such results may be 
of independent interest. 
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We start with a basic estimate for the case where /i = 0, n G 2N and the ceihng is 
the maximal possible configuration, A. 

Lemma 2.1. Let h = 0. For every a > there exists ci(a) > such that, uniformly 
in n £ 2N, k £ N and in the configuration a £ fl^^ of the floor, 

^A,-(yA(^) ^ i) ^ e-^^(»)^ (2.3) 

for every i £N. In particular, there exists p{a) > such that 

^A,a fU) ^^^yi^j^A ^ (2.4) 



Proof. By monotonicity, we can upper bound the probability in ()2.3p replacing a with 
V: 

^^-(y^(r?) ^i)^ ^^'^(y^(7?) ^ i). (2.5) 

Always by monotonicity, the right-hand side of (j2.5p is non-decreasing in n and k. 
Therefore, 

where M{v) is the number of plane partitions of volume v, and the right-hand side is 
just the limit n — )• oo, /c — )■ oo of the right-hand side of (|2.5p . The dependence on a, k 
and n has then disappeared. Since 

^W~°°^e.p(„,„^/») (2.7) 

for some (explicit) positive constants 01,02 llOj, one obtains immediately (j2.3p . Eq. 
()2.4p is obtained from ()2.3p just taking i = 1 and p(a) = 1 — exp(— ci(a)). □ 

Even if the ceiling ^ does not coincide with A, it is intuitive that the polymer rj^^^ 
gets more and more squeezed to ^ when k becomes large. This is one of the implications 
of the next result: 

Lemma 2.2. For every a > there exists 02(0) > such that the following holds. 
Uniformly in ^ ^ a £ in n,k £ N, \h\ ^ n and ^ x ^ n, 

/.-(t^O-) / ^^) ^ e-c2(a)(fc-i+i), (2.8) 
Moreover, for every < a < b < n, 

^^^^ (^ni^) / , Vo ^ X ^ ft) ^ g-C2(a)(b-a)_ (2.9) 

Proof. For a fixed 1 ^ x < n, let A G be defined as follows (see Fig. [3]): 

Ay := max(^^ - |y - x|, Vj^) , ^ y ^ n. (2.10) 

Note that ^ A and of course cr ^ V. Then, by monotonicity and the fact that = A^; 
we have 



«''^(#<e.)^^^'^(#<A.). (2.11) 
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Figure 3. Graphical construction of A. The thin full line denotes ^, 
the dashed line V and the thick full line is A, for a given value of x. The 
floor a is not drawn, since it has no influence on the construction of A. 

Now, rjx^ < Ax implies V^{ri) ^ k — j + 1. As in Eq. ()2.6p . from monotonicity it follows 
that 

and (I22D follows from {221) • 

Next, we prove ()2.9p . Denote by x^\ 1 ^ r ^ N{n, h) the positions of the particles 
of the ceiling ^, and let Xa^h '■= • • • ' ^i+mi those particle positions 

which are contained in the interval {a, . . . , 6 — 1} (the cardinality m + 1 of Ta^b does 

not exceed (6 — a) and can be zero). We use also the notation x^~^'^\ 1 ^ r ^ N{n, h) 
to denote the positions of the particles of the floor a. The event in the left-hand side 
of (j'2.9p implies that Xr'^^ > xl'^^ for every i ^ r ^ i + m. On the other hand, from (|2.8p 
we know that the event {x\''^^ = xf^^\/ 1 ^ j ^ A;} = {x^'^^^ = xf^^} has probability 
at least 1 — exp(— 02(0;)), uniformly in all the parameters. Assume that this event is 
not realized. In this case, the probability that {x^-'^^__-^^ = x^^j^_^ \/ 1 ^ j ^ k} = 
= xf^j^_-^} is again lower bounded by 1 — exp(— C2(a)): indeed, by monotonic- 
ity it is sufficient to consider the case where x^-;^^^ = xf^^ for every 1 ^ j ^ k, and to 



apply once more ()2.8p . Iterating this procedure, we see that the left-hand side of ()2.9p 
is lower-bounded by 

g-c2H|x„.,|_ ^2.13) 

One can then repeat the argument with the vacancies replacing the particles. The 
argument is the same except that vacancies have to be matched from left to right 
(while particles have been matched from right to left). Since the number of vacancies 
plus the number of particles in {a, ■ ■ ■ ,b} equals {b — a), one obtains immediately ()2.9p 
(modulo redefining 02(0)). 

□ 
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3. Proof of Theorem 11.11 

To prove the spectral gap estimate we shall use a three-fold decomposition that can be 
roughly described as follows. The first step, carried out in Section [3^21 allows to reduce 
the original process to a process with a bounded (i.e. independent of n, k) number s 
of polymers. This is achieved by estimating the spectral gap of a new process whose 
elementary moves consist in the updatings of sets of s full-length adjacent polymers. 
The second step, see Section 13.11 below, allows to further reduce the problem to a 
process where each of the s polymers has a bounded (i.e. independent of n, k) number 
r of particles. This is obtained by analyzing a coarse-grained dynamics where we update 
blocks of particles. The last step is a rough estimate for a system with s polymers and 
r particles, see Section [3.31 Finally, in Section [3.41 we prove the theorem by gathering 
all the pieces together. 

3.1. Particle block— dynamics. Let I'J^, for j < m, denote the equilibrium measure 
on E^^cr conditioned on the cr-algebra generated by the particle positions 

{x^^ , u = I, . . . , k} , V ^ j , and v ^ m . (3-1) 

As a convention, if v < 1 then we set (deterministically) Xv^ = for all j. Similarly, 
if V > N, we set x^^ = N. Recall that A'" = N{h, n) denotes the total number of 
particles. 

It will be convenient to have the following alternative notation. For a fixed integer 
£ we define the measures pf^ := fjl^-i j+^+i- Once the values of (|3.ip are given, pf^ 
is a distribution on the portion of paths rj^^^ in the segments {a^^^ -|- 1, . . . ,b^^^ — 1}, 
j = 1, . . . ,k, where a^-'^ := x[''}^_-^, b^^^ := x^l^^^-^^, see Figure IH 




Figure 4. An illustration of the definition of the measure pfg for the 
paths from Figure[TJ Here ^ = 2, i = 5. The shaded region is the portion 
of paths distributed according to p^^, while the rest of the configuration 
is frozen. 

We will study the following Markov chain with state space -E^.^, for a given pair 
^, o" E ^^in) with ^ ^ a. This auxiliary process, which we call the particle block- 
dynamics, is reversible w.r.t. the equilibrium measure ^^'^ . 

We have N{h, n) independent Poisson clocks with parameter 1. When the i~th clock 
rings we consider the current configuration ry, and update the portion of paths rj^^'^ in 
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the segments {a^^^ + 1, . . . , b^^^ — 1}, j = 1, . . . ,k with a sample from the conditional 
distribution pf^. The rest of the configuration rj is left unchanged. In other words, we 

(i) (i) 

are removing from the system all particles at positions xy_^, . . . , x^^^, j = 1, . . . ,k and 
we are replacing them by a sample from pf^^. 

The generator of this process can be written as 

N 

Of = J2[pfAf^-f] ' (3-2) 

i=l 

where / denotes a function / : E^^a- — ^ ^ and pfiif) is the function rj ^ J fiOpfiidC) 
for X = x{ri). Since pf^ are conditional expectations we see that p^''^ (^fipfeif) ~ /)) ~ 
p^''^ (^{pf^if))'^ — Pfeif'^i^i so that the Dirichlet form of the process is 

N 

- P^-" ifGf) = t"^'" (Va^Pf,. (/)) • (3-3) 

i=l 

In particular, Q is self-adjoint in L'^{p^''^). Let gap(^) denote the spectral gap of this 
process, which of course depends on the choice of ^ G N. 

Proposition 3.1. For any a > and any A; G N, there exists i = i{a, k) such that 
uniformly in h, n, and £,,a ^ ^i,n '^^ have 

gap(a) ^ 1 . 

To prove Proposition 13.11 we use a coupling argument. Consider two evolutions 
rj{t),ri' (t), t ^ 0, of the Markov chain described above, with initial conditions rj and rj' 
respectively, where rj,rj' £ E^^a- A well-known argument (see e.g. Proposition 3 in [9]) 
shows that, for any coupling P of the two evolutions 



gap(^) ^ — lim inf - log 

t— >-oo t 



max 



X P (7?(t) / rj'it)) 



(3.4) 



Therefore, to prove Proposition 13.11 it is sufficient to establish Lemma 13.21 below. 



Lemma 3.2. For every a > and A; S N there exist £ G N,7 > and a coupling P 
of {ri{t),ri' (t)) such that, uniformly in the starting configurations and uniformly in the 
parameters h,n and ^,0" G 

F{r]{t) / r]'{t)) ^ kN{h, n)e^" e"* , t ^ . 

Proof. Consider the distance 

N k 

i=i j=i 

(i) (i) 

where x- , y^- denote the positions of the i-th particle of the j-th polymer in the 
configurations and rj' respectively, and 7 > is to be determined later. Note that the 
minimal non-zero value of d^(-, •) is e~^'^. Therefore, by Markov's inequality 

nvit) ^ v'{t)) ^ e'y^Eid,{^{t),l^'m . 



11 



We need to show that we can define a couphng such that 

E{d^{T]{t),r]'{t))) ^ kN{h,n)e-* . 

From path couphng ([I]), see Lemma [331 below for the details, it will be sufficient to 
exhibit a Markovian coupling such that 

^ - d^iv,v'), (3.6) 

t=o+ 

for all pairs r],rj' satisfying do{T],r]') = 1, i.e. when there is only one discrepancy in 
the particle configurations. In this case, d.y{rj,rj') = e~^^ where i is the label of the 
discrepancy. 

Next, we specify the Markovian coupling, and prove that it satisfies ()3.6p . If we use 
the same Poisson clocks for the two evolutions, the infinitesimal generator Q of the 
coupled dynamics can be written as 

N 

Q = ^(i^;,,, - 1) , (3.7) 

i=l 

where -Ej^ denotes a (not yet specified) coupling of the local equilibria pf^ of the block 
of 2^ + 1 particles around the i-th particle for the k polymers, cf. (|3.2p . 

To prove ()3.6p . we may choose the coupling in such a way that if 77,77' have a 
single discrepancy at a given particle label % and at a given polymer label j, then 
[-E't,^£d^](r7, 77') = for all v such that v — i^i^v^l. By construction, there are at 
least I such blocks for any fixed % (this is the case if e.g. i = 1 or i = A^). 

Let E~ = £'j_£_i^^ (respectively, = £'j+£+i^^) denote the coupling corresponding 
to the block of 2£ + 1 particles just to the left (resp. just to the right) of particle i. Note 
that if e.g. i ^ i + 1 then there is no block just to the left of i and we may set = 1 
for such i. Similarly, ii i ^ N — i then there is no block just to the right of i and we 
can set E^ = 1 in this case. Since all other blocks give a trivial contribution to (|3.6p 
we see that 

[§d^]{v,v') (3.8) 
^ - (77, 77') + [{Er - l)d^](7?,7?') + [{E+ - l)d^] (77,77') . 

Recall that d^{r],r]') = e~^^. Then we can estimate 

k i+2e+l 

[E^d,]{r,,v')=Yl E '~'''^t 

u=l v=i+l 

00 _ 

^ A: e-^^ ^ e-^^ = d,{v, v') , (3.9) 

1 1 c 

11=1 

where we have bounded by 1 the probability of a discrepancy. 
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On the other hand, denoting by Fy the event that there exists u £ {1, . . . , A;} such 



that / y\'^^, we have 



i-l 

1, 

t 

u=l v=i-2e-l 

21+1 21+1 

v=l v=l 

with the convention that Ip^ = i — v < 1. 

Below, see discussion after ()3.1ip . we prove that for a suitable choice of the coupling 
there exists c = c(a, k) > independent of 7 such that that 

^^[l^J^e"'^ l^vi^2i + l. (3.10) 

If we assume this estimate, from (|3.8p . (j3.9|) we conclude that if e.g. 7 = c/2 then 
(j3.6p follows for £ sufficiently large (depending on k and a). In particular, the proof of 
Lemma 13.21 will be completed once we prove (|3.1Up . 

We turn to the proof of ()3.10p . Let x, y denote the collections {x^^}, {yi^^} of all 
positions of particles of two configurations 77, r]' G E^^^- and consider the associated 
probability measures i^fmi^^m defined in (|3.ip . For a fixed pair of integers j, m, let 
^ _ p,m j-jg independent coupling of j and That is, we are 

freezing all particles labeled i^j — m — 1 or i^j and we sample the m particles 
labeled i = j — m,...,j — 1 according to the independent coupling of j and 

^J-m-i j- ^' y agree up to j — m — 1 if = yi^\ for all u = 1, . . . , A; and 

for all V ^ j — m — 1. 

We claim that there exists e > depending only on a and k such that, if x, y agree 
up to j — m — 1, then 



u X 



yY\ , for all u = l,...,k) ^ e . (3.11) 



Let us first show that (j3.1ip implies (|3.10p . First of all, let us sample x)^"i\,yl^\^, 
u = 1, . . . ,k, using with j = i and (supposing for simplicity i > 2i -\- 1) m = 21 + 1 
(recall that i is the index appearing in the proof of Lemma l3.2p . Then (|3.1ip implies 

that there is a full matching x'f^-^^ = y'f^i for all n = 1, . . . , A; with probability at least e. 
Thus ()3.10p holds when v = 1 and e~'^ = 1 — e. The case > 1 is obtained by recursion. 
Namely, the coupling E~ can be further defined as follows. If we have a full matching 

x[^i = y^i for all u = 1, . . . , A; then we can match all the remaining particles. If we 
do not have the full matching we sample particles labeled i — 2 hy V with j = i — 1 
and m = 2£ and a suitable choice of the values x[^i,yl^^ (the ones that were sampled 
in the first step). Again, if we have a full matching for particles labeled i — 2 we can 
match all the remaining particles labeled i — 3, . . . ,i — m. We repeat this procedure at 
later steps. Since ()3.1ip shows that there is a probability at least 1 — e~'^ to have a full 
matching at every step the bound ()3.10p follows. 
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To prove (j3.1ip we are going to use the same argument as in the proof of Lemma 

see the proof of (j2.13p in particular. We first sample the pairs xi^^ , yi^^ , r = 
j — m, . . . , j — 1 with an independent coupling for the corresponding marginals. Note 
that even if x^p 7^ yj^^ the positions x^-^^, have the same ground state (i.e. minimal 
position) dictated by the ceiling ^ and the (common) boundary conditions for the 
particles labeled 1, . . . , j — m — 1 of polymer 1. From the proof of Lemma 12.21 we know 
that this implies that the matching event x^^]^^ = Vj^li has probability at least 6"^ where 

6 = 6{a) > is the probability that x^j'^-^ equals its ground state position. Next, we 
sample x^\yi'^\ r = j — m, . . . , j — 1 with an independent coupling of the marginals 
conditioned on the configuration of {xi^\yi^^}r=j-m,...,j-i which we extracted in the 
previous step. We claim that, conditionally on the occurrence of the matching x^^\ = 

Vj^li, the matching event x^-^^ = yf^i has probability at least 5^. Indeed, even if the 

(2) (2) 

ceilings "felt" by the two copies of the positions xf- ,yr of particles of polymer 2 are 

(2) (2) 

in general distinct, the ground state positions for x-_^,y-J_^ are dictated only by the 

positions which are now assumed to coincide. In particular the argument 

from the proof of Lemma 12.21 again applies. This procedure can be repeated until the 
last polymer is reached and in conclusion the independent coupling gives probability 
at least e = 5'^^ to the full matching event in p. lip . This ends the proof of Lemma 

□ 



In the proof of Lemma 13.21 (cf. ()3.6p ) we have used a continuous time version of the 
so-called path coupling argument, see [1] or e.g. [HI Theorem 14.5] for the usual discrete 
time version. For the sake of completeness we give a proof in the next lemma. 

Lemma 3.3. Suppose that 13. 6\) holds for all pairs rj^rj' such that dQ{rj,rj') = 1. Then 
K{d.y{r]{t),ri' (t))) ^ kN{h,n)e^* , for all t and for all initial data rj^rj' . 

Proof. We may define a graph having as vertices the elements of E^^^j by declaring a 
pair r], to be an edge whenever c?o(??, ??') = 1. For any ry, ry', let ^^(r/, r/') denote the set 
of paths connecting r/ and ry', i.e. w G 0,(ri,ri') if a; = (wi, . . . ,0Or), ooi = rj,uJr = rj' and 
do(wj, Wj+i) = 1, i = 1, . . . , r — 1. Then one checks that 

r-l 

d-fiv^v') = min ^^7(^4, a;i+i) . (3-12) 

u)en{ri,'r)') ^ 

To prove (j3.12p observe that it suffices to exhibit one path which achieves the equality 
since by construction it is clear that dy{rj,rj') ^ l^[=i d^i^i,^i+i) for any u S Q{r],r]'). 
Such a path can be informally defined as follows. Consider the bottom paths ry(^), ry'*^'^) 
and the positions x where r/i'^'' > rjx''^ . Then move one by one the particles of tj^^^ in this 
region until we have rj^^ ^ rj'x^'^ everywhere. Then consider the paths 7^^^^^), ry'^'^"^) 
and the positions x where r]^ > 77^^ ^\ As before, move one by one the particles 
of r]^^~^^ in this region until we have r]i^ ^ rjx^ everywhere. We repeat this 
procedure until we reach the top paths. At this point we have reached a configuration 
rj such that rj ^ r] (everywhere). Next we start from the top paths rj^^\ri'^^^ and 
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consider the region where ryx < v'x ■ We can move one by one the particles in this 
region until we have ry*-^^ = rj'^^^ (everywhere). We repeat with the paths r/^^^ry'^^^ and 
so on until we reach the bottom polymers labeled k. This construction produces a path 
which realizes the minimum ()3.12p since we never used more than the strictly necessary 
moves. 

From the triangle inequality, for each pair of initial conditions rj, rj' we have 

r-l 

E[d^{l^{t),7]'m ^ ^E[d^(Ci(t),Ci+i(i)], 

1=1 

where we call Ci , . . . , Cr the minimizing path in (|3.12|) and Ci (t) , . . . , Cr (t) the corre- 
sponding trajectory. In particular, subtracting d-y{r],rj') = YliZi d-y{Ci,Ci+i)j dividing 
by t and letting t ^ we obtain 

r-l 
i=l 

Since each term in the r.h.s. above is of the form ()3.6p with do{Ci,Ci+i) = 1 we have 
that the assumptions imply 

[gd^]{ri,ri') ^ - d.,{ri,ri') , 

for arbitrary initial conditions. Therefore (f'{t) := E[d^(?7(t), r/'(t))] satisfies ^^{t) ^ — 
(p{t), which implies the claim since (p{0) < max^^^i d^{ri,ri') < kN{h,n). □ 

3.2. A dynamics with fuIl-poIymer moves. The next ingredient which enters the 
proof of Theorem 11.11 is a dynamics where each move consists in updating (2s + 1) 
whole polymers, s G N. As usual, we let G (with a ^ and /i^''^(-) denotes 
the law on E^^^ for k polymers with floor and ceiling a,(,, as defined in (jl.3p . To each 
1 ^ i ^ ^ is assigned an independent Poisson clock of mean 1. When the clock labeled 
j rings, we update the polymers 7?^"^ with index max(j — s, 1) ^ u ^ min(j + s,k), 
sampling the new configuration according to the law 

Z>^.,(.) := fj,^'^ I ^(max(j-s-l,0))^^(min(j+s+l,fc+l))^ ^ (^g^^g) 

with the convention that r/^*^) := ^ and := a. Call the generator of this 

dynamics. 

Proposition 3.4. For every a > there exists s := s{a) G N such that, uniformly in 
n,k, \h\ ^ n and on o" ^ ^, one has 

gap(A^) ^ 1. (3.14) 

Proof. The general structure of the proof is similar to that of Proposition 13.11 but the 
coupling argument is rather different. Given p > 0, we define the distance function 
Dp{-, •) by setting for every rj,rj', 

DM) :=E^"^'^E 2 • ^'-''^ 

i=l a:=l 
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In analogy with (j3.12p . one checks that 

r-l 

Dp{r],r]')= min y^Dp{uJi,uJi^i), (3.16) 

where in this case one requires that L'o('^i; f^j+i) = 1 for i < r. 

Given two initial conditions (f/, f/'); ^^t {r]{t),r]' (t)) be the corresponding evolutions. 
As in the proof of Proposition 13. H cf. Lemma 13.31 to prove ()3.14p . it is then sufficient 
to prove that for every a > there exists a choice of s G N, /) > and a coupling P of 
{r]{t),ri' (t)) such that 



^-D,{r],r]') (3.17) 



whenever DQ{rj,rj') = 1. 

Let {fj, f)') satisfy the latter condition, with the single discrepancy consisting in f/i*^ = 
f)x'^^ +2. Li analogy with (|3.7p . we write the generator of the coupled dynamics as 

k 

M:=Y, {Ej,s - l) , (3.18) 
i=i 

where Ej^g is a coupling (to be specified) of Oj^s for the two configurations. For all j 
such that i — s^j^i + swe can choose the coupling such that [Ej^gDp] (v, v') = 0. 
Moreover, if j G {1, . . . , k} \ {i — s — 1, . . . , i + s + 1} we can choose the coupling such 
that [Ej^sDp]if],v') = Dp{f],f]'). One has therefore 

[MDp]if],fj') (3.19) 

< - sDp{f],fi') + [(^.+.+1,. - l)Dp]{fi,fi') + [(^.-.-1,, - l)Dp]{fi,f]'). 

It is clear that the last two terms may be non- negative, and that they vanish if i ^ k — s 

or i ^ s + 1, respectively. 

Let us analyze first the easier case of Ei^s+i,s in which case, it is worth recalling, we 

are updating polymers labeled i + l,...,i + 2s + l. Since f/^*) ^ f]'^^^ while i^(«+2'5+2) _ 

^/(i+2s+2)^ by monotonicity there exists a coupling Ei^g^i^s such that one has rj^^^ ^ ij'^^^ 

for every i<j<i + 2s + 2. Moreover, since r/(*) and f/'^*) differ only at x, we can 

(7) '(7) 

choose Ei^s+i,s such that rjy = rjy for every i<j<i + 2s + 2 and y outside the 
interval {a_, . . . , a+}, where a+ := inf{y > x : rjy^'^'^^'^^^ = f]y^} and a_ := supj?/ < x : 
^(i+2s+i) _ j^W|_ consequence, going back to the definition of Dp{-, •), 

00 

[Ei+s+i,sDp]{fj,fi') ^ e-^'i'i+s+hs ((a+ - «-)') • (3.20) 

j=i+i 

Thanks to ()2.9p there exists 03(0) > such that for every u > 

i)i+s+iAa+ - a_ = n) ^ g'^^^")". (3.21) 
From this one deduces immediately that there exists 04(0, p) < 00 such that 

[Ei+,+i,sDp]{fi,fj') ^ C4{a,p)Dp{n,f)'). (3.22) 
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Finally we deal with Ei-s-i,s- We have from (|3.15p 



[E, 



i—s—l,s 



i—s—l,s 



'n-1 



iJ) _ 



y 



(3.23) 



j=i-2s-l 

Again, we can choose the coupling such that tj'^^'^ ^ ry^-'^ and rjy^ = rjy^^ for y ^ 

{&_,..., where 6_ = sup{y < x : rjy'' ^'^ = fjy "^^ ^^}, and similarly for 6+. In 
analogy with ()3.2ip one has 

Then, 

n-l 



y=\ 



Using & and (ICTD 

n— 1 



i— s— l,s 



n— 1 



i— s— 1 



5_)21 g-C4(Q)(i-j) 




Therefore, going back to ()3.23p . 

oo 

[E,.s-i,sD,]{v, f/O ^ c,{a)D,if], fj') e-('=^(")-'')^ = c,{a)D^^^) (f/, ,7'), (3.26) 



r=l 



where we chose p := /3(a) = C4(a)/2. From ()3.19p . ()3.22p and ()3.26p one concludes that 

[MDp^^)]ifj,f]') ^ - is- C7(a))I)p(«)(r/,f/'). (3.27) 

At this point, it is sufficient to choose s := s{a) := [07(0)] + 1 to get (|3.17p . □ 

3.3. An estimate for k polymers with r particles. The last ingredient we need 
for the proof of Theorem 11.11 is a rough estimate on the spectral gap for a system with 
k polymers, each with r particles; it is important that this bound is independent of 
the lengths ni,. . . ,nk of each polymer. Consider a configuration rj E E^^^j C fl^^ with 
k polymers, each with N particles. Let x = x(rj) denote the corresponding particle 
configuration. Fix 0^j<m^N + l and consider the probability measures vj^ 



defined in ()3.ip . If we freeze all particles labeled i ^ j and i ^ m we can perform 
the local-update dynamics defined in ()1.9p for the r := m — j — 1 particles labeled 
J + 1, . . . , m — 1. This process is clearly reversible with respect to fj^^- Dirichlet 
form is given by 



^f,mif) = Yl E ^j:™ [Var„,.(/)] 



U = l (u) , 

x=x^ +1 



where Var„^^(/) has the same meaning as in (jl.Op . As usual, below we use the notation 
Var^(/) = i^(/^) — i^(/)^ for any probability measure u. 
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Lemma 3.5. For every a > 0, € N and r G N, there exists c = c(a, fc, r) suc/i t/iat 
/or any pair j, m with m — j — 1 = r, for any choice of all other parameters and for 
any f : J^^_„ ^ M 

q,m{f)-^c\^i,^jf). 

Proof. The only delicate point here is that the length of the portion of paths where the 
r particles live is arbitrarily long and we need an estimate which does not depend on 
that. We introduce a further family of measures as follows. Let ^^'^ denote the law of 

r/(") according to the measure vj„^ conditioned on the value of the paths 77^*^ , i ^ u. If 
we let 

x("'-l 
x=xf + l 

denote the corresponding Dirichlet form we know from the k = 1 version of Theorem 
11.11 (see the remark following the statement of Theorem II. ip that 

^^";::(/)^'^(«) Var,x,.(/), (3.28) 

for some constant 5{a) depending only on a. Taking expectation w.r.t. i^Jrn^ using 
f^^[f^^((7)] = vjrnid) function g, and summing over u in ()3.28p . by definition 

of Sf^mif) ^^'^^ 

k 

£l^{f) ^ 5{a)Y,^l 

u=l 

Therefore it remains to prove that for some c = c(a, r, /c) > one has 

k 

^cVar,x^(/), (3.29) 

n=l 

for all functions /. To prove this estimate we observe that the left-hand side of (|3.29p 
coincides with the Dirichlet form of the Markov chain described as follows. Attach 
independent rate 1 Poisson clocks to the labels u = l,...,/c. When a label u rings 
update the whole path between x^f^ + 1 and x^m — 1 according to the distribution vf'^ 
(that is, freeze all other polymers and update polymer u with a sample from i^j'm)- 

Thus, the following rough coupling argument will suffice for the proof of ()3.29p . 
Namely, consider the Markov chain started in the minimal configuration (i.e. each of the 

k polymers starts in the minimal path compatible with the particles x^"^ and Xm^). Let 
Et denote the event that up to time t the Markov chain has never visited the maximal 
allowed configuration. It is not hard to prove a bound of the form P(i^'f) ^ c~^e~'^^ for 
some c = c(a, r, k) > 0. This, in turn, implies the desired spectral gap estimate ()3.29p 
using monotonicity and the bound ()3.4p . 

To prove F{Et) ^ c~^e~'^^ we may reason as follows. Consider the event Fg that in 
the time interval [s, s + 1) the following sequence of k successive updatings appears: for 
i = 1, . . . ,k the i-th update is for polymer u = i and the configuration sampled from 
lyj'^ is the maximal allowed path for the i-th polymer given the current configuration. 
Since there are r particles, an application of Lemma 12.21 shows that P(-Fs) ^ p where 



Var^x,« (/) 



Var^x,„(/) 
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p = p{a, r,k) > depends neither the lengths = Xm — Xj of the polymers, nor 
on the configuration of r] at time s. By construction, on the event Fs we know that 
the Markov chain visited the maximal allowed configuration at least once in the time 
interval [s, s + 1). Therefore the event Et implies that none of the events Fs occurred 
for s = 0, . . . , [t - IJ . This implies F{Et) ^ (1 - pf'^ ^ (1 - p)"^e-P*. The proof of 
Lemma 13.51 is complete. □ 



3.4. Putting everything together: proof of Theorem ll.il Once Proposition [3TT1 
Proposition 13.41 and Lemma 13.51 are established, the proof of Theorem 11.11 is obtained 
through a chain of comparison inequalities. Indeed, Proposition 13.41 mav be restated as 



Var(/)^ ^f^[y^^r,^M)] > 



(3.30) 



where Var(/) is the variance w.r.t. fi := fi^''^ and uj^s is the conditional probability 
measure defined in (|3.13|) . Here s = s{a) is a fixed parameter (2s + 1 is the number 
of polymers to be updated at each step in the process with generator A4 appearing 
in Proposition 13. 4t they are the polymers labeled u, with max(j — s — 1, 1) < u < 
min(j + s + 1, k)). For each j in (|3.3U|) we apply the estimate of Proposition 13.11 with 
the number of polymers equal to k = s. This yields 



N 



Var;>.^(/)^ ^i>j,s Varpx^(/) 



(3.31) 



i=l 



where pf^ is the conditional probability measure appearing in (j3.3|) and £ = £{a) is a 
fixed parameter {2£ + 1 is the size of the block of particles to be updated at each step 
in the process with generator Q appearing in Proposition 13. ip . Finally, recalling the 



definition pf^ 



^ (cf. ()3.ip ) and applying the bound in Lemma 13.51 (with 



r = 2£ + 1 and A; = s) we know that for every i in (|3.3ip we have the estimate for all 
X = x(r/): 



-1 



j + S + l -^i^l + l' 

Varpx^(/)^c-i E pfA^^^nAf)] 

U=j-S-1 ^_„(") ,1 



(3.32) 



Note that the constant c = c(a, s, £) depends only on a here. Also, note that since 
there are at most 2£ + 1 blocks of particles covering a given position x in the path we 
obtain 

N j+s+1 n-1 

J2 [Varpx^(/)] ^ c-i(2^ + 1) E E ^^'S [Var„,x(/)] , (3-33) 

u=j—s—l x=l 



i=l 



where we use the property of conditional expectation Uj^s pfiid) 
any function g. 



'^j,sig) valid for 
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Putting everything together and using fi [^'j,s(g)] = n{g) for any g, from (j3.30p - (j3.33p 
we see that 

k j+s+1 n— 1 

Var(/)^c-i(2^+l)^ ^ ^ ^ [Var„,,.(/)] 

j=l u=j—s—l x=l 

k n— 1 

< c-i(2^ + l)(2s + 1) ^ ^ [Vari,,(/)] 

i=l x=l 

= c-\2i + l)i2s + l)£if). 

This shows that gap(a) ^ c'(q), with c'(a) = c{2i + l)~^(2s + which completes 
the proof of Theorem ll.li □ 

4. Proof of Theorem 11.21 

For clarity of exposition we give the proof first of all in the simpler case where the 
ceiling ^ is the maximal configuration, A, in Later, in Section [4.31 we sketch the 

modifications which are needed in the general situation. 

4.1. The case of maximal ceiling. Fix some positive T > 0. By monotonicity and 
the definition (jl.lip of mixing time, it is clear that a sufficient condition for T^[x ^ T is 
that the first (random) time when the dynamics started from the minimal configuration, 
jjyO) = = I, . . . ,k}, hits the maximal configuration, {rj^^^ = = 1, . . . ,k}, is 
smaller than T with probability at least 1 — l/(2e). Again by monotonicity, it is easy 
to convince oneself that this random time does not decrease if one replaces h with zero, 
n with 2M := 2max(n, k), k with M and a with V. 

Therefore, to prove Theorem 11.21 (in the case of maximal ceiling) it is sufficient to 
prove the fohowing. Let h = 0, n = 2M, A: = M, ^ = A, cr = V, cf. (fLT]), dLl]); start 
the dynamics from the minimal configuration := {rj^j^ = V for all j = 1, . . . , Af}, 
and call t{M) the first time when the maximal configuration := {r]^^^ = A for all 
j = 1, . . . , M} is reached. Then, 

Theorem 4.1. There exists C{a) > such that with probability larger than 1 — l/(2e) 

t{M) ^ C(a)M (log M)^ (4.1) 

Theorem 14.11 will be proven just after the statement of Theorem 14.21 below. Before 
stating Theorem 14.21 we need to introduce in Sections 14.1.11 and 14.1.21 a few auxiliary 
definitions. 

A notational convention When in the following we say that an event occurs "with 
large probability" (w.l.p.), we mean to say that the probability of the complementary 
event goes to zero for M — )• cxd faster than any inverse power of M. Since, as will 
be clear, we have to exclude only polynomially many (in M) events which w.l.p. do 
not occur, by the union bound we have that the occurrence of at least one of these 
events still goes to zero for M — )• oo. For simplicity of exposition, and when there is no 
risk of confusion, we will often just pretend that an event which occurs w.l.p., occurs 
deterministically. 



20 



PIETRO CAPUTO, FABIO MARTINELLI, AND FABIO LUCIO TONINELLI 



4.1.1. From polymer configurations to subsets of the cube. The proof of Theorem 14.11 
becomes more intuitive if one interprets a configuration rj £ 2m ^ ^ subset s of 
the cube Cj\/ := [0, Af]^ C M'^. This mapping is just another way to see the mapping, 
introduced in Section 11.1. 11 between polymer configurations rj and boxed plane parti- 
tions. To define precisely this correspondence, divide first of all Cm into elementary 
cubes of unit side, which we label with the (integer) coordinates r = (ri, r2, rs) of their 
point of smallest norm (observe that ^ < M). Then, to a given rj G ^m2M 
we associate s = s(ry) C Cm, a union of elementary cubes, by establishing that the 
elementary cube labeled r belongs to s if and only if (cf. Figure [5]) 

v'iit?L<M-n-r2. (4.2) 




Figure 5. A graphical construction of the set s{r]). Given r] G ^%2M^ 
in order to obtain s{r]) do the following, for rs = 0, . . . , (M— 1): a) draw 
the configuration of jy(*^~''3) (thick line in the left drawing); b) rotate 
the picture by Stt/A anti-clockwise and shrink it by a scale factor 1/^/2. 
The shaded region coincides with the horizontal section of s{r]) at height 
r^ < h ^ r^ + 1. It is obvious from this construction that the subset 
s{rj) thus obtained is a monotone subset of the cube Cm- 



It is easy to check that: 

• if rj ^ rj' then s' C s (note that the inequality is reversed!) 

• the maximal configuration r]~^ defined above corresponds to the empty subset 
of the cube: s~ = 0. Conversely, the minimal polymer configuration rj~ corre- 
sponds to the maximal subset s"*" = Cm, i-e., the full cube. We will sometimes 
refer to s~ as to the ground state, for obvious reasons. 

• for every r/ G ^^m2M' '*(^) ^ monotone subset of Cm, i-e., if the elementary 
cube labeled r belongs to s, then so do also all the elementary cubes r' such 
that r[ ^ ri,r2 ^ r2,r'^ ^ r^. 

The equilibrium measure is described in terms of s by: 

p-2q|s| 
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where |,s| denotes the number of elementary cubes contained in s, i.e., its volume. 

The dynamics of Section [1.21 can also be explicitly described in terms of s. Here, let 
us simply remark that the elementary moves of the Markov Chain consist in adding 
or removing a single elementary cube, with the constraint that s remains a monotone 
subset of Cj\/ after the update. Observe also that, under the ry -f-)- s correspondence, 
the upward drift felt by the polymers rj^^^ during the dynamics translates into the fact 
that the upper boundary of s feels a drift in the direction (—1, —1, —1). 

4.1.2. An auxiliary dynamics. To avoid a plethora of [-J , we assume that (logM)2 £ 2N 
and that K := (M/(logM)^) G N. Divide Cm into sub-cubes (called blocks from 
now on) of side (logM)^, indexed by f = {vi,V2,V3) with ^ Vi < K , and such that 
the point of B^ with minimal norm is (logM)^f (of course, Vi are integers). Given 
V, we will also define B^,B~ to be the half-blocks obtained cutting By horizontally 
into two equal parts {B^ will denote the bottom one). Call st the configuration at time 
t, which starts from the completely full configuration s'^ at time t = 0. 

The idea behind Theorem 14.11 is to consider an easier dynamics st such that st C st 
almost surely and, calling t(M) the first time st reaches the empty configuration s~ , 
to show that t{M) satisfies (|4.ip with probability at least 1 — l/(2e). The statement of 
Theorem 14. II then follows immediately by monotonicity, since t(M) < t{M). 

Let 

r := r(M) := C(a)(log Mf/7, (4.4) 

where C{a) is the same as in (|4.1|) . We define now two deterministic sets C Cm 
which (roughly speaking) coincide with at t = 0, are empty after time C(a)M(log M)^, 
and such that st satisfies C st C , w.l.p. and for all t ^ (the latter property 
is non-trivial and it is the content of Theorem 14.21 below). 

{S^}t ^ 0) is defined as follows (see also Figure [6] for a graphical definition): 

• Sq = Cm 

• if i G N and (i — l)r < t ^ ir, then contains all and only the blocks By such 
that V satisfies 

V3 + 2{vi+V2) i^5{K -l)-i. (4.5) 
On the other hand, letting for ease of notation 

V,^:=CM\St, 

{Sf^}t > is defined as follows (see also the caption of Figure [6]): 

• S^cS+cCm 

• li By G K-[t/Tj By-.(;- G ^-[t/Tj Icast one choice of i = 1,2,3, then 
By G (cj are the canonical base vectors of Z^, and [x\ := max{n G Z : 
n ^ x}) 

• U By G ^t^j/t-j but there is no i = 1,2,3 such that By^ei £ V^^^/^j , then 
B+ G y+ but B- G St- 

The following properties of are immediately checked: 

• is the union of blocks By, while is the union of blocks and of half-blocks 

By 
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• Sr is left-continuous, non-increasing, constant in the time intervals ((i — l)r, ir] 

• is right-continuous, non-increasing and constant on intervals [{i — l)r, it) 

• are monotone subsets of Cm for every t. 

• For t ^ (6/7)C(a)M(log = 6Kt, is empty and contains only B~ ^ 

{K := M/(logM)^ and r = C(a)(log M)®/7 were defined at the beginning of 
this section). 

As an example, if < t ^ r then = Cm \ B(^k-i,k-i,k-i)i while if 2r < t ^ 3r 
then 

Sf = Cm \ (-B(js:-2,x-i,ii"-i) U B(^k_i x-2,k-i) uf=i -B(j£--i,K-i,i<:-i))- 



• z 




, X 



Figure 6. The set for 14r < t ^ 15r. Small cubes denote blocks By 
of side (logM)^, and there are K = M/(logM)^(= 10) of them along 
each side. At time intervals of 2r, a new diagonal set of columns with 
vi + V2 = const starts to move downwards: from then on, it moves one 
block down each time interval r. Once a column is empty, it stays empty 
forever. Roughly speaking, the set contains all the blocks which are 
below a plane perpendicular to the vector (2, 2, 1) and which moves at 
constant speed, of order (logM)~^, in the direction (—2, —2, —1). The 
set S^' for ir ^ t < {i + l)r can be obtained simply by taking 
for some {i — l)r < t ^ it and adding a half-blocks B~ on top of each 
incomplete but not empty column, and also to each empty column which 
is adjacent to a non-empty one {Sf~ is not drawn in the picture). 



We define the auxiliary dynamics st by establishing that it has the same law as st 
conditioned on the event that S~ C s„ for every u ^ . Remark that by monotonicity 
we can couple st and st in such a way that st C for every t ^ 0; also, remark that 
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once the deterministic set is empty there is no more constraint on the dynamics 
(which does not mean that st = st after that time!). 

The basic point is the following: 

Theorem 4.2. W.l.p., for every ^ t ^ one has 

Sr CstC 5+. (4.6) 

Of course, the lower bound is trivial by the very definition of st (and holds not only 
w.l.p. but with probability one). As a side remark, it will be apparent below that it 
would be sufficent to have the above statement with replaced by M^'^'' for some 
e > 0. 

Proof of Theorem \4.1\ (assuming Theorem \4.2\) . Thanks to ()4.6p . for all times 

(6/7)C(a)Af(logM)'' ^ t ^ 

one has st C = B^^^^y which is a subset of a cube of side (logM)^. This is 

just because S^^ = for t ^ (6/7)C(a)M(log M)^, as we observed just after the 

definition of . By point (1) in Lemma 14.31 below, this implies that, within time 
(6/7)C(a)M(logM)6 + 0((logM)6) ^ C{a)M{logMf < M^, st has hit the ground 
state s~ at least once, with probability at least 1 — l/(2e). Since st C st, this implies 
that St has also hit s~ and ()4.ip follows. □ 

4.2. Controlling the auxiliary dynamics Sf. The first ingredient of the proof of 
Theorem 14.21 is the following lemma, which gives a rough upper bound on the mixing 
time in a cube of size M: 

Lemma 4.3. For every a > there exists Ci(a) < oo such that for every M G N 

T^i.i^Ci{a)M^. (4.7) 

Moreover, there exists 6*2(0) > such that for every T > the following is true: 

(1) the probability that t(M) ^ TM^ is smaller than exp{—C2{a)T) . 

(2) with probability at least 1 — T e~^^^°'^^^ , the volume of st is at most M/10 for 
all times t{M) ^t^T. 

Proof. We know that 

where J\f{v) was defined after formula (j2.6p and is just the number of plane partitions of 
volume V. Then, it follows from ()1.13p (and the fact that the gap is uniformly positive) 
that the mixing time is 0{M^). 

Prom ()4.7p it is immediate to deduce (modulo redefining Ci(a)) that the probability 
that t{M) ^ Ci{a)M^ is greater than some e(a) > 0. (Indeed, Lemma 12.11 implies 
that there exists a set A of configurations such that /i^'^(A) > 1/2 and such that 
all the configurations s £ A can be reached via at most m = m{a) Markov Chain 
moves from the ground state s~ , for some m(a) < 00 independent of M). From this, 
one easily deduces that, for every n G N, the probability that t{M) > nCi{a)M^ is 
smaller than (1 — e(a))", i.e., claim (1). Indeed, if the evolution has not hit s~ before 
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time (n — 1) Ci{a)M^, just restart the dynamics from the maximal configuration at 
t = in — l)Ci{a)M^: this can only make t{M) larger, by monotonicity. 

To prove statement (2), observe first of all that for all times t > t{M) (or, more 
precisely, conditionally on t{M) < t) the distribution of st is stochastically dominated 
by the equilibrium distribution ^jl^'^ . On the other hand. Lemma l2. II tells us that 

/i^'^(|s| ^ M/10) ^ exp(-ci(a)M/10). (4.9) 

Secondly, the number of Markov Chain moves in the interval [0, T] is a Poisson random 
variable C, with average TM"^ , and an elementary computation shows that for a Poisson 
random variable C,\ of parameter A one has 

P(Ca ^ ra) ^ e-"('°s("/^)-^). (4.10) 

We have therefore, calling ti,i = 1, . . . , ^ the random times when the updates occur, 

¥{3t G [t{M),T] : \st\ > M/10) ^ e"^"^ (4.11) 

+P (C ^ AM^T; 3i : t{M) ^T,\st^\-^ M/IO) 

^ e"*^'^ + 4M2 T exp(-ci(a)M/10), 

where in the last inequality we used the union bound and (j4.9p . □ 

Notational convention: in the rest of this section, we will use for simplicity of exposi- 
tion expressions like "for all times larger than to" to mean "for all times to ^ t ^ M^". 

For the next lemma we need some notations. Let ^ j < 2{K — 1) and let 5(j) 
be the set of configurations s such that C s if {vi + V2) < j and i?^, n s = if 
(fi + V2) > j + 1 or if f 1 + t;2 G + 1} and ^ 2. Let s^^^ be the maximal 
configuration in 5(j), see Figure [7| (of course, both S{j) and s^^^^ depend on M). We 

define a dynamics {s[''^}t ^ by requiring that (in law) it equals our usual dynamics, 

with initial condition s^IIq = s™^^ and conditioned on the event that G '^(j) for 
every t ^ 0. 

Lemma 4.4. W.l.p., the following holds for all times t ^ r. If {V1+V2) = j and = 1, 
or if {vi + V2) = i + 1, then n s^^ = 0. // {vi + V2) = j and V3 = 0, B+ D s^^ = 0. 

Proof. As in the proof of Theorem 14.11 we introduce an auxiliary dynamics {s[''^}t ^ 
for which the claim is easier to prove, and such that almost surely s[''^ C s[''^ (so that the 
claim follows also for s^^). The auxiliary dynamics is defined simply by requiring that 
its law equals that of {s^^jt > conditioned on the event that, for all times ^ t ^ r/2, 
B-^ C sp^ if (^1+^2) = j and vi E 2N + 1. In other words, such blocks B^ are frozen and 
remain completely full up to time r/2. This implies that during this time interval the 
blocks which are not frozen evolve independently by groups of at most six, see Figure [HI 
Applying Lemma 14.31 one sees that at time r/2 the situation is the following (w.l.p.): 

(1) By n 4"/2 = if (^^1 + V2) = J + 1 or {vi + V2) = J>i e 2N and = 1; 
the same remains true for all later times (of course here we are using the fact 
that, since r = C(a) (log M)^/?, the probability that t{{logM)^) ^ r/2 is 
0(exp(— C"(a)(log M)^)), cf. statement (1) of Lemma l4.31 i.e., w.l.p. such event 
does not occur.) 
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Figure 7. The maximal configuration s^^x _. S{j). The dark 

region is the one which is constrained to remain completely full during 
the evolution of sf , while the dashed one remains empty. The white 
region is the one which can evolve. Note that, for graphical convenience, 
the axes are drawn with orientations which differ from those of Fig. [6l 



(2) If (ui + V2) = j, vi G 2N and = 0, the volume of By n s^J^ is at most 
(logM)^/10 (in particular, n s^^j^ = 0); the same remains true for all later 

times. Here we are using statement (2) of Lemma l4.3t with T = M^. 

(3) all blocks marked by a star in Figure [8] are of course still completely full. 

During the interval (t/2,t] the blocks marked by a * in Figure [8] are free to move. 
However, now the pairs of blocks vi + V2 = j,vi £ 2N + l,t>3 E {0,1} evolve inde- 
pendently for vi different, thanks to point (1) above. Therefore, applying once more 

Lemma 14.31 and point (2) above, one sees that for t ^ r one has By n = for 

vi+V2= j, f 1 G 2N + 1, = 1, and B+ D sp^ = for vi + V2 = j, vi £ 2N + 1,V3 = 0. 
The claim is proven. □ 

Proof of Theorem \4.2[ It is sufficient to prove the following: for every i G N U {0} such 
that ir ^ M^, one has (w.l.p.) 

St C for every t ^ ir. (4.12) 

For i = the statement is trivial since Sq = Cm- 

Now we assume that the claim is true up to a certain i, and we show that it holds also 
for i + 1. It is convenient to introduce the definition of column Cy^^y^{H) with (integer) 
base coordinates Q ^ vi,V2 < K and height ^ H ^ K: this is just a parallelepiped 
of height (log M^H whose base is the square of side (log M)^ such that the base point 
of minimal norm has coordinates ((log M)^f 1, (log M)^f2). Note that, for each t, 
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Figure 8. Here we drew only the region where is ahowed to evolve 

(the white region in Fig. [7]). Under the dynamics s[''\ blocks marked by 
a -k are constrained to remain full up to time r/2. As a consequence, up 
to that time the remaining blocks (shown in gray) evolve independently 
by groups of at most six. 



can be viewed as composed of such columns: in the case of 5^" the heights 
H^{vi,V2) take values in {0, 1, ... , K} while in the case of Sf' we call them H^{vi, V2) 
and they take values in {0, 1/2, 1, . . . , K — 1/2, K}. Note also that, by construction, 
H^{vi,V2) depends on (fi,f2) only through vi + V2 (see Figure E]). 
To complete the inductive proof, we need to prove that for all {vi, V2) 

[st n C^,,V2{K)] C C^^,v2{H^i+i)r(.vi,V2)) for every t ^ (i + l)r. (4.13) 

The following cases can occur (keep Figure [6] in mind) : 

(1) H^.^-^^^^{vi,V2) = K. In this case, (j4.13p is obvious because also H^^^-^^^{v\, V2) = 

(2) 1 < H^^^^^_^{vi^V2) < K. In this case, we have that 

Cv„v,{H[,^i)A^uV2)) C [st n C,,,v,{K)] C C,uv2{H^i+i)M,V2) + 2) (4.14) 

for every ir < t < (i + l)r. The lower bound is trivial by the definition of 
the dynamics st, while the upper bound follows from the inductive hypothesis 
()4.12p and from the fact that the definition of implies that ^ H^{vi,V2) — 
Hf (vi,V2) < 2. Since we want to prove (j4.13p . by monotonicity we can assume 
that all the columns labeled {wi,W2) with wi + W2 < vi + V2 are completely 
full during the time interval ir < t < (i + l)r. But then, as we shall argue in 
a moment, in the time interval (ir, (i + l)r) the column st n C^^^^2(K) evolves 
independently of all the others, and an application of Lemma [4.31 implies ()4.13p . 
since H'^_^_-^^{vi,V2) is just H7._^-^.{vi,V2) + 1/2, cf. the caption of Figure El 
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To see that the column (^1,^2) evolves independently of all the others in the 
interval ir < t < (i + l)r, note that it can be influenced only by the columns 
labeled {wi,W2) with W1+W2 = f 1 + ^2 + 1- However (cf. Figure [6|) in this case 
-f^(7+i)r(^i' "^2) = 1^2) — 2 so that, from the induction hypothesis (cf. 

{vi,V2)). (4.15) 

In other words, the column {'Wi,W2) is too low to influence the column (fi,V2)- 

(3) H^._^-^y{vi,V2) = 1. Again one has ()4.14p and one can assume by monotonicity 
that the columns with W1 + W2 < vi + V2 are completely full in the time interval 
ir < t < (i + l)r. The argument proceeds like in the previous case once one 
realizes that for W1 + W2 = V1 + V2 + I one has H^{wi,W2) = 1/2, so that by the 
induction hypothesis Sj fl C^tj^^^jl-^) C«ii,w)2(l/2) for t> ir and such column 
cannot influence the one labeled (^1,^2)- 

(4) It remains to consider the case of the columns with H^^-^^-^^ = 0. Define 

j := max {wi + W2 : HT.^. {wi,W2) > 0} + 1 < 2i^ - 1, 

with the convention that j := if the set is empty. It is convenient to distinguish 
two sub-cases: 

(a) if wi + ^2 > j + 1 then, by definition of S^^ one sees that H^{vi,V2) = 0, 
so that (|4.13p follows (both sets are empty) from the induction hypothesis 

dm. 

(b) if j < V1 + V2 < J + 1 then by monotonicity we can assume that all columns 
with W1 + W2 < j are completely full in the time interval it < t < (i + l)r, 
and on the other hand we know that all columns with W1+W2 > j + 1 are 
completely empty for t > ir. Also, we know from (|4.12p that, always for 
t ^ ir, [st n Cwi,w2i^)] C Cwi,w2i'^) for {1^1,102) such that wi + 11)2 G 

+ 1}. We can therefore apply Lemma 14.41 to deduce that for all times 
t ^ (i + l)r the columns with wi + W2 = j + 1 are completely empty, 
while if wi + W2 = j then [st H Cwi,w2iK)] C B~{'Wi,W2,0). Recalling the 
definition of (in particular, as explained in the caption of Figure E]) we 
have therefore proven (|4.13p for all columns {'Wi,W2) such that wi + W2 € 

□ 

4.3. The general case. Here we prove Theorem 11.21 in the general situation where 
^ is not the maximal configuration in The proof is conceptually similar to the 

one where the ceiling is maximal, and therefore some arguments will be only sketched 
(see however Remark 14.61 below, where an important difference between the two cases 
is pointed out). 

As in section WA\ we can assume by monotonicity that h = 0, n = 2M, k = M and 
o" = V. First of all, it is important to realize that the maximal configuration (ground 
state) {t]^^^ = ^ for all j = 1, . . . , M} corresponds to the subset of Cm defined by the 
property that the elementary cube labeled r = (ri, r2, ra), with G {0, . . . , M — 1} (cf. 
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Section 14.1. ip belongs to s if and only if 

^M-n+ra < M - ri - r2 (4.16) 

(we still call the ground state s~, even if it is no longer the empty set as in previous 
section). We note also that the equilibrium measure (|1.4p is given in this case by 

„-2a\s\s~ 

^^''"(^)= J2,^^_,-2aW\s-\ ''s'^s-. (4.17) 

Of course, the sum in the denominator is only over the configurations s' <Z Cm which 
are given by unions of elementary cubes. The dynamics st just coincides with that 
described in Section I4.1.H conditionally on the event that st D s~ for every t (of 
course, the initial condition st=o has to verify the same property). 

As in Section we cut Cm into K-^ blocks of side (logM)^ G N. Consider 

the dynamics started from some st=o D s~ . The basic estimate which allows to prove 
Theorem 11.21 is the following: 

Proposition 4.5. For every a > there exists C{a) < oo (independent of M, ^ and 
St=o) such that the following holds w.l.p.: 

stC{s-UAi) (4.18) 

for every {6/7)C {a)M (log ^ f ^ where Ai = A2LI A3, 

A2 ■= D{By : V3 = 0, By D s^ contains at least one elementary cube} (4.19) 

and 

A3 := U{By ■.V3 = 0, 5(^i_i,„2,o) ^ (Cm \ ^2) and ^(^^^^^.^o) ^ (Cm \ ^2)}, (4.20) 
see Figure\M 

Note that, in the case ^ = A, one has A2 = 0, A3 = -B(o.o,o) ^-nd Proposition 14.51 
follows from Theorem 14.21 above. 

Proposition 14.51 is proven below, and now we show that it does imply Theorem 11.21 

Proof of Theorem \1.2\ (for general ceiling ^ ) assuming Proposition \4.5l 
Note that the volume of Ai satisfies 

\Ai \ s-\ ^ AM(logM)^, (4.21) 

so that 

min (s\s- CsC{s"U Ai)) ^ e-c8(«)M(iogM)4 ^^^22) 

s:s-CsC(s-UAi) 

for some cg > (we used the fact that the number of configurations s satisfying 
s~ C s C (s~ U ^i) is smaller than 

2l^iV-l : 

this would be the exact number of configuration if there were no monotonicity con- 
straints on s). Call, for ease of notation, 

Tm := C(a)M(logM)^ (4.23) 
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Figure 9. The sets Ai and s~ seen from above. Squares should be 
imagined to have side (logM)^, and here K = M/(logM)^ = 8. The 
dark region is a horizontal section of s^, while the lightly colored (re- 
spectively, the dashed) region is the set A2 \ s~ (resp. ^3) seen from 
above. s~ extends vertically up to height M, while ^2)^3 extend only 
up to height (log M)^ (i.e., one single block). Note that the block i?(o,6,o) 
belongs to A^^ because neither 5 q) nor -6(0,5,0) belong to Cm \ A2 
(in particular, B(^_i q q-j falls out of Cj\/). A similar remark holds for 
-6(7,0,0)- 

From (jl.lSp , ()4.22p and Theorem 11.11 one easily deduces that the mixing time of the 
dynamics constrained to s~ C st C [s~ U Ai] is 0(M(logM)^). Since {6/7)Tm + 
M(logM)^ < Tm < M2, Theorem O easily follows. □ 

Remark 4.6. It is important to notice a crucial difference between the proof of The- 
orem [L2] for maximal ceiling A and for generic ceiling ^. In the former case, we have 
shown a stronger statement, i.e. (cf. Theorem 14. ip that the dynamics started from 
an arbitrary st=o hits the ground state s~ within a time of order M(logM)^. This is 
related to the fact that the ground state s~ has an equilibrium weight which is pos- 
itive, uniformly in the system size. In the generic case, i.e. when the ceiling ^ has a 
jagged shape, there is no reason to believe that the dynamics hits s~ within such time; 
indeed, its equilibrium weight (an therefore the inverse of its hitting time) may well be 
exponentially small in M. 

4.3.1. Proof of Proposition By monotonicity, it is clear that it suffices to prove the 
claim for st=o = Cj\/ . As in the proof of Theorem 14. H we introduce two deterministic 
subsets of Cm- If are the sets which were defined in Section 14.1. 2[ then we 
establish that 



St ■■= St Lis , 



(4.24) 



30 



PIETRO CAPUTO, FABIO MARTINELLI, AND FABIO LUCIO TONINELLI 



while 

S+:=5+UAi. (4.25) 

We note that from the discussion of the properties of in Section |37L2] it follows that 
for t ^ {6/7)Tm one has = s~ and = Ai. The claim of the proposition then 
follows if we can prove, in analogy with Theorem I4.2t that w.l.p. one has 

St C 5+ (4.26) 

for every ^ t ^ M^. The proof of this fact proceeds with the help of an auxiliary 
dynamics st, which dominates st, and whose law is that of st conditioned on the event 
that {su D for every u ^ 0}. Since the proof is very similar to that of Theorem 14.21 
we do not give details. The only fact which requires some care is that we cannot apply 
Lemma 14.31 to the blocks which have a non-empty intersection with s~ , since a certain 
number of its elementary cubes are frozen to be full for all times. 

The extra result we need concerns therefore the evolution of a single cube of side 
(logAf)^ and with arbitrary ceiling. Let | G ^?,2(iogM)2' ^~ ^ C(iogM)^ = ^(o,o,o) 
be the ground state corresponding to the ceiling ^ and st be the evolution started from 
the full configuration i5(o,o,o) and constrained to 

S~ C St C 5(0,0,0) 

for every t > 0, and call Tmix its mixing time. Needless to say, its invariant measure 
is /i^'^. Notice that when s" = i3(o,o,o) when = V) the dynamics is trivial 

{st = -8(0,0,0) fo'^ all times) while when s~ = (i.e., when = A) the forthcoming 
lemma is already implied by Lemma 14.31 (iust replace M with (logM)^ there). 

Lemma 4.7. For every a > there exist C4{a) < oo and C^i^a) > such that 
uniformly in M and ^ 

fmi.^Ci{a){\ogMf. (4.27) 
Moreover, with t defined as in (j4.4p . 

P (3 r ^ t ^ : (st \ S-) n S+ q / 0) ^ ^-C,ia){logM)\ ^^ ^S) 

What Eq. (j4.28p is saying is essentially that for all times larger than (log M)'^Tjnix 
but smaller than the upper half of the cube under consideration contains only the 
elementary cubes which are imposed by the constraint s~ G st- Useless to say, this 
means that w.l.p. the event in the left-hand side of (|4.28|) does not occur. 

Using Lemma [4. 71 as a substitute for Lemma [4. 3 1 the proof of Proposition 14.51 is easily 
concluded in analogy with the proof of Theorem 14.21 We spare the reader additional 
details. □ 

Proof of Lemma \4. 71 The proof of (j4.27p is completely analogous to that of (j4.7p and 
uses the fact that the spectral gap of the dynamics is positive, uniformly in M and ^. 
To prove ()4.28p . let us recall the well known inequality which relates the total variation 
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distance from equilibrium of a reversible Markov Chain at time t with its mixing time 

mix • 



sup||P/(-)-M- 



(4.29) 



where is the state space of the Markov Chain, fi its invariant measure and Pf the 
law at time t, if the initial condition at time zero is s. Call ti,l ^ i ^ the random 
times when Markov Chain updates occur in the time interval [r, M^], and observe that 
in our case ^ is a Poisson random variable of parameter (M^ — r)(logM)^. One has 
then, using (|4.1U|) . that the left-hand side of (j4.28p is upper bounded by 



C ^ 4M2(log M)4 and 3 i ^ C : {h, \ 5") n ^,0) ^ 



(4.30) 



On the other hand, calling Pt{-) the law of st and defining the set 
U:={sC 5(0,0,0) : (« \ s') n S+ ^ / 0}, 

one has 



ist\s )ni3+o,o)^ 



Pt{U) ^ Pt{U)-fi^'''{U) +/z«>(f/) (4.31) 



In the last step, we used (j4.29p for the first term and Lemma 12.21 to estimate the 
equilibrium probability of U. The claim ()4.28p then follows from (|4.30p . a union bound 
and the fact that r/Tmix ^ C6(Q;)(log M)^ for some CQ{a) > 0. □ 
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